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ABSTRACT 


The  determination  of  the  critical  speeds  of  a  flexible  rotor  is 
usually  carried  out  on  the  assumptions  of  negligible  bearing  mass  and  stiffness. 
In  this  report  the  effect  of  bearing  mass  and  stiffness  on  the  natural  frequencies 
of  a  uniform  flexible  shaft  is  found.  Two  specific  cases  are  considered:  the 
shaft  of  symmetric  cross-section  rotating  in  bearings  of  different  horizontal 
and  vertical  stiffness,  and  the  unsymmetric  shaft  rotating  in  symmetrically 
stiff  bearings.  The  frequency  equations  for  various  types  of  end  constraints 
are  solved  by  means  of  a  simple  graphical  procedure.  It  is  found  that  bearing 
mass  and  stiffness  usually  have  considerable  influence  on  all  but  the  first  or 
second  critical  speeds. 
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NOMENCLATURE 


Sj.S^  principal  flexural  rigidities  of  shaft  (S2  S^) 

S^S^+S^/Z  mean  flexural  rigidity 

L  length  of  shaft 

mass  per  unit  length  of  shaft 
_Q_  angular  speed 

X  natural  frequency  (radians  per  second) 

k  flexibility  of  bearing 

M  mass  of  bearing 

displacement  of  shaft  in  first  principal  direction 

V2  displacement  of  shaft  in  second  principal  direction 

Z  distance  along  shaft 

Non-Dimensional  Variables  and  Parameters 
T  =  t  time 

u  =  Vj/L  displacement 

v  =  V^/L  displacement 

z  =  Z/L  variable  distance  along  shaft 

£  =  (S2-S^)/(S2+Sj)  measure  of  shaft  anisotropy 

p  ■  nV  'h  L^/S  angular  speed 

h  =  y  gL  /S  reduced  gravity  force 
3 

K  =  kL  /Sratio  of  bearing  and  shaft  flexibilities 

m  =  M/  y  L  ratio  of  mass  of  bearing  to  mass  of  shaft 
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DYNAMICS  OF  FLEXIBLE  ROTORS 


OBJECTIVES 

Because  of  its  importance  in  rotating  machinery,  the  dynamics  of 
flexible  shafting  has  been  the  object  of  much  theoretical  and  experimental 
study.  This  research  has  not  only  attempted  to  determine  the  vibrational 
behavior  of  various  systems  as  a  function  of  operating  speeds,  shaft  para¬ 
meters  and  support  conditions,  but  also  devised  methods  of  minimizing  these 
motions.  For  a  review  of  previous  analysis  and  balancing  techniques,  the 
reader  is  referred  to  Armour  Research  Foundation  Report  No.  K206  and  the 
extensive  bibliography  contained  therein. 

Despite  the  most  exacting  balancing  procedures,  there  usually 
exists  a  small  deviation  of  the  center  of  gravity  of  the  shaft  from  the  axis 
of  rotation.  As  the  shaft  rotates,  this  mass  eccentricity  produces  a  peri¬ 
odic  centrifugal  force;  it  is  by  this  excitation  that  flexible  shafts  are  brought 
into  unwanted  vibration.  Generally,  these  oscillations  are  negligible  but  at 
certain  speeds  of  rotation,  the  system  may  vibrate  at  a  level  which  is  intol¬ 
erable  for  proper  operation  of  the  machinery  or  its  environment.  In  addi¬ 
tion  to  causing  structural  damage,  excessive  vibration  may  interfere  with 
delicate  instruments,  be  of  discomfort  to  humans,  or  generate  noise.  The 
problem  of  noise  reduction  is  of  particular  importance  in  naval  construction 
since  sound  detectors  are  a  primary  means  of  locating  ships  above  and  below 
the  surface. 

The  initial  step  in  developing  intelligent  means  of  vibration  attenua¬ 
tion  is  an  analysis  of  the  motion.  Once  an  understanding  of  the  interrelation 
of  various  system  parameters  has  been  made,  the  designer  can  rationally 
proceed  to  a  system  of  desired  smoothness. 

In  standard  analyses  of  rotating  shafts,  it  is  assumed  that  the  bear¬ 
ing  mass  itself  has  negligible  effect  on  the  characteristics  of  the  vibrations. 
This  is  certainly  the  case  in  a  rigid  bearing;  in  practice,  however,  this 
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idealization  is  sometimes  impossible  to  achieve  because  of  practical  limi¬ 
tations  on  the  size  of  the  supporting  structure.  Indeed,  it  is  possible  that 
a  bearing  is  rigid  in  one  direction  and  quite  flexible  in  the  other  direction. 

The  Bureau  of  Ships,  Department  of  the  Navy,  Washington  25,  D.  C.  , 
established  this  research  project  under  Contract  No.  NObs -86805  in  order 
to  investigate  the  effects  of  heavy  flexible  bearings  on  the  dynamic  behavior 
of  flexible  rotors.  It  was  believed  that  this  information  would  be  of  value 
to  anyone  concerned  with  rotating  machinery  as  a  means  to  understand  and 
avoid  unwanted  vibrations. 
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INTRODUCTION 


We  shall  investigate  the  dynamics  of  a  flexible  shaft  rotating  at 
constant  speed  in  heavy  flexible  bearings.  To  keep  the  problem  as  general 
as  possible  without  obscuring  it  with  mathematical  complexities  we  shall 
assume  that  the  shaft  has  different  stiffnesses  in  two  mutually  perpendicular 
directions  but  is  otherwise  uniform  along  its  length.  The  bearing  masses 
and  flexibilities  may  all  be  different. 

Several  authors  have  considered  the  dynamics  of  non-symmetric 
n 

rotors'  ;  all  of  these,  however,  have  not  included  the  effect  of  bearing 
mass  and/or  flexibility.  Foote,  Poritsky,  and  Slade^  have  treated  a  sys¬ 
tem  comprised  of  a  single  mass  on  a  weightless  non-symmetric  shaft  sym¬ 
metrically  mounted  in  massless  flexible  bearings.  Their  investigation  in¬ 
volved  two  ordinary  coupled  differential  equations  with  periodic  coefficients. 
The  results  indicate  that  shaft  anisotropy  alone  will  cause  instability  over  a 
range  of  speeds  while  shaft  and  bearing  anisotropy  can  cause  unstable  oper¬ 
ation  over  as  many  as  three  ranges  of  speed.  Furthermore,  they  noted 
that  the  inclusion  of  damping  does  not  necessarily  remove  these  instabilities. 
Because  of  the  nature  of  the  differential  equations,  the  mathematics  is  quite 
lengthy  and  it  is  doubtful  whether  the  approach  could  be  modified  to  include 
the  effect  of  massive  bearings.  Moreover,  since  they  considered  only  a 
one-lump  parameter  system,  there  can  be  no  information  as  to  stability  of 
the  continuous  system  at  super-critical  speeds. 

Kellenberger^  has  treated  the  uniform  non- circular  shaft  rotating 
in  immovable  bearings.  The  partial  differential  equations  of  motion,  when 
referred  to  a  rotating  coordinate  system,  are  easily  solved.  His  method, 
however,  is  applicable  only  to  systems  whose  boundary  conditions  are  in¬ 
variant  under  rotation.  It  is  found  that  there  occur  an  infinity  of  critical 
speed  ranges  that  vary  with  the  amount  of  anisotropy  of  the  shaft.  Second¬ 
ary  critical  speeds  induced  by  gravity,  are  also  determined. 


*  Numbers  in  brackets  denote  references  collected  at  end  of  this  report. 
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In  this  report,  the  vibration  of  an  undamped  shaft  in  heavy  self- 
aligning  flexible  bearings  is  studied.  A  method  of  characteristic  exponents 
is  used  to  investigate  the  stability  of  the  system.  Expressions  are  found  in 
closed  form  for  the  dividing  surfaces  between  stable  and  unstable  domains. 
The  secondary  critical  speeds  are  found  to  bear  a  simple  relation  with  one 
of  these  expressions.  The  method,  being  sufficiently  general,  is  applicable 
to  other  types  of  end  constraints. 

As  a  special  case  of  the  above  analysis  the  frequency  equations  for 
various  practical  round  shaft  systems  are  found  and  the  critical  speeds  are 
presented  in  graphical  form  as  functions  of  bearing-shaft  mass  and  stiffness 
ratios. 
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PROBLEM  STATEMENT 


Let  the  shaft  of  Fig.  1  have  unequal  principal  stiffnesses  S ^  and 
S^.  As  the  shaft  rotates,  its  stiffness  in  any  fixed  coordinate  system  will 
be  a  periodic  function  of  time  completing  two  cycles  for  every  revolution  of 
the  shaft.  Consequently,  the  equations  of  motion  in  any  fixed  coordinate 
system  will  have  periodic  time -dependent  coefficients.  When  referred  to 
rotating  coordinates  fixed  in  the  shaft  along  the  principal  directions  of  stiff¬ 
ness  all  coefficients  in  the  equations  of  motion  are  constant.  If  the  hori¬ 
zontal  and  vertical  flexibilities  at  the  bearings  are  unequal,  the  flexibilities 
in  the  rotating  system  are  also  periodic  functions  of  time.  In  order  to  avoid 
periodic  coefficients  in  the  boundary  conditions,  we  shall  assume  that  the 
bearings  possess  equal  flexibility  in  two  mutually  perpendicular  directions. 
This  assumption  facilitates  the  analysis  of  the  non-symmetric  shaft  and  is 
unnecessary  in  the  treatment  of  the  isotropically  stiff  shaft.  We  are  thus 
able  to  investigate  two  cases;  the  anisotropic  shaft  in  isotropic  bearings  and 
the  isotropic  shaft  in  anisotropic  bearings. 

So  that  the  results  may  be  applicable  to  any  system,  the  analysis  is 
carried  out  non -dimensionally.  In  terms  of  u  and  v,  the  dimensionless 
displacements  in  the  rotating  system,  and  other  dimensionless  variables  and 
system  parameters  defined  elsewhere,  the  equations  of  motion  for  the  aniso¬ 
tropic  shaft  are 


,,  .  ^u  .  2  .  d^u  ,3v  . 

‘  1  57  p  -u-2rr>  ‘ 


-  h  sin  t  +  p  f  (z)  (1) 


^4  P  *2  ^  . 

0+0  — z  +  P  (  -  v+2~)=-h  cost  +  p^  g{z)  (2) 

dz*  dr*  3 


Here  the  parameter  «  £l)  is  a  measure  of  the  shaft  asymmetry  with 

vanishing  £  implying  an  isotropic  shaft.  The  angular  speed  of  the  shaft 
is  characterized  by  p  while  h  measures  the  gravity  force.  The  functions 
f(z)  and  g(z)  represent  the  mass  eccentricity  of  the  shaft. 

We  shall  investigate  the  motion  of  a  shaft  mounted  in  heavy,  self- 
aligning  bearings.  These  lead  to  boundary  conditions  commonly  referred 
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to  as  pinned -pinned  ends;  the  analysis  that  follows  could  apply  equally  as 
well  to  any  other  type  of  end  constraint.  The  equations  of  motion  of  the 
bearings  are 


■*3 

(1-.)  |4 

#z 

+  Kj  u  t  mjp2 

>2 
/  d  u 

— *3?> 

2 

+  Kj  v  +  m^p 

-~2 

-v  +2  p 

-.3 

-(I-.)  L» 

dz5 

+  K2  u  +  m2p2 

.  52U 

-u  -2  - — 

♦<u«>4 

dz3 

+  K2  v  +  m2p2 

.  0  v 
3tZ 

.  ->  Du 

"v+2  TF 

(3) 


The  additional  conditions  for  pinned-pinned  end  constraints  are 

2  2 

at  z  =  0,  1  ^-4  =  0  =  0  (4) 

3  a 

Equations  (1)  -  (3)  are  derived  in  Appendix  A. 

The  problem  is  to  find  the  displacements  as  functions  of  the  bear¬ 
ing  masses,  m^,  m^,  stiffnesses  Kj,  K2>  speed  p,  and  asymmetry  factor 
€  .  For  certain  values  of  speed  and  asymmetry  we  should  expect  a  steady- 
state  solution  while  at  other  speeds  the  displacements  will  increase  without 
limit.  While  infinite  amplitudes  will  not  occur  physically,  the  system  will 
operate  with  excessive  vibrations  at  these  so-called  critical  speeds.  The 
first  task,  then,  is  to  determine  these  critical  speeds. 
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THE  HOMOGENEOUS  SOLUTION 


Because  of  the  nature  of  the  non- symmetric  shaft,  the  homogeneous 
solution  is  of  considerable  importance.  Even  though  damping  may  be  pres¬ 
ent  in  the  shaft,  the  homogeneous  part  of  the  solution,  while  being  transient, 
can  still  become  unbounded.  This  seemingly  paradoxical  situation  is  ex¬ 
plained  in  Appendix  B. 

Since  the  equations  of  motion  (1),  (2)  have  constant  coefficients,  the 
homogeneous  solution  can  be  expressed  as 


u  =  A  e 


ocz  +  i  p  t 


B  e 


oc  z  +  i  p 


(5) 


where  A,  B,  oC  and  p  are,  in  general,  complex  constants  to  be  determined 
by  the  end  conditions.  When  equations  (5)  are  substituted  into  (1)  and  (2),  it 
is  found  that  these  constants  are  related  by 


(l  -«)<*4  -  Pz  (1  +  |i2) 


p  (1  +  p  )|  A-2p(i  iB  =  0 


2p  \x  i  A  + 


-  “ 

{  l+c)(X4  -  p2  (  1  +  |12)  B  =  0 


(6) 


For  a  non-trivial  solution,  the  determinant  of  coefficients  of  (6)  must  vanish. 
Performing  this  operation,  we  obtain  a  relation  between  oi  and  p  which  is 

4 

quadratic  in  oc  .  Since  the  time  exponent  p  determines  whether  the  mo¬ 
tion  is  stable,  it  is  advisable  to  eliminate  OC  in  favor  of  p  .  Thus 


oc  =  |p2/(i-«2)J  jji+i*2)  +  V‘W>2+  4^2<1-‘2)  J 

For  each  possible  p  there  are  eight  values  of  OC  and  sixteen 
constants  of  integration  Aj,  Bj  which,  from  (6)  and  (7)  are  related  by 


(7) 


Bi  -  -i0i  aj 


(j  =  1,  2,  3,  . . .  8) 


(8) 


where 


Oj  =  -«(i+p2)+yr 


2  ,,  .  2.2  ,  .  2 2. 

(1+p  )  +  4p  (1-*  ) 


J^p(l+c)j 
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The  upper  sign  is  to  be  taken  for  j  =  1,  2,  3,  4  and  the  lower  sign  for 
j  s  5,  6,  7,  8.  The  solution  (5)  may  now  be  written  as 


oo 


*  ■  E  .  E  ^ 

root  p  j  =  1 


OC  ^(p)  z  +  ipT 


(9) 


00 


__  O'  j(lO  z  +  i|*T 

v=  -i  2_  OjV 

root  p  j  =  1 

Here  the  summation  is  extended  over  all  characteristic  values  p  which 
shall  be  determined  from  the  boundary  conditions.  For  each  as  yet  un¬ 
known  p  the  constants  of  integration  are  found  from  the  boundary 
conditions  (3)  and  (4),  which  in  terms  of  the  notation  (8)  and  (9)  become 

(1-c  )  *.3  +  Kj  +  mlP2  (-p  2  -  1  +  2i0.)l  =  0 


j  =  i  L 

)  Aj  -(l+«)  -  KjiO.  +  m^2  (iO^p2  +  i»j  +  2ip  )j 

'j>J  =0 


=  0 


8 

E 

j  =  i 
8 

E 

j  =  i 


(10) 


Ot 


A.e 

J 


-(l-« )  o-3  +  K2  +  m2p2  (-p  2  -  1  +  2i0 


(1+0  iOj°j ^  -  K2  iOj  +  m2p2  (iOjP2  +  iOj  +  2ip)  =  0 


1 


and 


y~  oc.z  a.  =  o 

L —  J  J 

j  =  1 

8  OL  j 

2  J 


w 

Ee.  oi}  a.  =  o 
j  j  j 

j  =  i 


E°v 


j  =  i 


e  A.  =  0 
J 


8  OL  . 

2  J 


y~  o.  oi, 
L —  j  j 
j  =  i 


A.  =  0 
J 


(H) 


ARMOUR  RESEARCH  FOUNDATION  OF  IlllNOIS  INSTITUTE  OF  TECHNOLOGY 


-  9  - 


ARF  Final  Report 
No.  K274  . 


The  frequency  equation  for  p  is  obtained  by  setting  the  determinant 

of  coefficients  of  (10)  and  (11)  equal  to  zero.  In  this  manner,  a  non-trivial 

solution  for  A.  will  be  obtained.  For  any  system  configuration,  all  para- 
J 

meters  except  p  ,  of  course,  are  known.  When  all  frequencies  are  real, 
the  vibrations  are  steady-state;  when  the  frequencies  are  complex,  the  mo¬ 
tion  is  transient.  The  transient  will  decay  or  become  unbounded  depending 
upon  the  imaginary  part  of  the  complex  frequency. 

It  is  evident  from  Eq.  (9)  that  the  system  is  stable  for  given  system 
parameters  p,  t ,  m,  and  K  whenever  real  (ip  )  is  non-positive.  It  is  a 
simple  matter  to  show  that  the  frequency  equation  is  an  even  function  of  p 
so  that  both  p  *  and  -p  *  are  roots.  If  p  *  is  a  complex  root  and  real 
(ip*)  is  non-positive  then  real  (-ip*)  is  positive  so  that  the  amplitudes 
may  increase  exponentially  with  time.  Clearly,  stability  is  assured  when¬ 
ever  all  frequencies  are  real.  Since  instability  occurs  whenever  the  fre¬ 
quencies  become  complex,  the  bounding  surfaces  between  stable  and  unstable 
domains  are  the  loci  of  parameters  p,  c  ,  m,  and  K  for  which  p  =  0  is  a 
root  of  the  frequency  equation.  Physically  p  =  0  implies  there  is  no  time 
variation  of  the  displacements  in  the  rotating  coordinate  system.  The  shaft, 
then,  is  performing  a  "frozen  whirl"  at  the  onset  of  instability. 

As  indicated  in  Appendix  C  when  p  approaches  zero,  the  eighth 
order  determinant  reduces  to  the  product  of  two  determinants  of  fourth 
order.  The  bounding  surfaces  are  found  from 
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The  loci  of  parameters  for  which  Aj  or  A2  vanishes  represents  the  be¬ 
ginning  or  end,  respectively, of  an  unstable  domain  (as  the  speed  p  is  in¬ 
creased).  As  «  approaches  zero,  these  two  curves  in  the  p-t  plane 
converge  to  the  natural  frequencies  of  vibration  of  the  symmetric  shaft. 

The  determinants  (12)  can  be  expressed  as 

A,=(l-€)2ar,6(l-cos£^cosh<xri)+Q!ri3(i-£)  [(K1+K2)-(m1+m2)( l-t  )4^*|  x 

x  jsinh^jcos d^-coshA^sinA^”|  +2  Jm ^(l-c  |^m2(  1-c  sin^sinh^sO 

A2=(l4-€  )2^56(l-cos^5cosh«^)+^53(l4«  )  |(K1  +  K2)-(m1+m2)(l+€  )<*5*J  x 

x  J^sinhd^gcosodf^- coshes  in  +  2  |m  ^(l+t  Jm2(  l+«  )^^-K^Jsin°^sinh«^=0 

For  the  special  case  of  symmetric  mounts,  m1=m2=m,  Kj=K2=K 
the  loci  reduce  to 

( 1-*  )2  ^^(l-cosO^cosho^H  2(  1  -t  )Ad/13(K-mp2)(cos<^19inhod'1-sin<^1cosh«^1) 


+  2  (K-mp2)2  sinA^sinh^j  =  0  (13a) 

and 

(1+t  )2c^5°(l-cosO'5cosh«'5)+  2(  l  +  «  Jod'g3 (K-mp2) (cos  ^sinhA^-sin^cosh^) 


t  IV 

w  «A*  ••***  / 


o  i  nVi  o/ . 


(13b) 


The  unstable  speed  ranges  for  a  pinned-pinned  shaft  in  immovable  bearings, 
found  previously  by  Kellenberger,  are  obtained  from  Eqs.  (13)  when  K 
approaches  infinity.  This  results  in 


sin  d^sinhq^  =  0 

sin^t-sinhQ^g  =  0 
Q^,g  —  nr 

or  _ 

(nr)2  yi  -  «  <  p  <  (nr)2  yTTT  (n=l,2 . )  (14) 

ARMOUR  RESEARCH  FOUNDATION  OF  ILLINOIS  INSTITUTE  OF  TECHNOLOGY 


-  13  - 


ARF  Final  Report 
No.  K274 


EFFECT  OF  GRAVITY  -  SECONDARY  CRITICAL  SPEEDS 


The  gravity  solution  can  be  treated  completely  and  independently 
of  the  previous  solution  for  the  homogeneous  case.  Let  us  assume  a  par¬ 
ticular  solution  of  Eqs.  (1)  and  (2)  in  the  form 

u  =  x  (z)  sin  t  v  =  y  (z)  cos  t  (15) 

With  this  formulation,  the  time  dependence  can  be  factored  out  of  the  equa¬ 
tions  of  motion  leaving  the  following  equations  in  the  amplitudes  x(z)  and 

y  (®): 


(1-0  — r-  +  2  p2  (-x  +  y)  =  -h 

dz* 

,4  , 

(l+€ )  — £  +2pZ(x-y)  =  -h 

dz 

Eliminating  y  (z)  from  Eqs.  (16)  we  obtain 


(16) 


y  (*)  = 

-  f(l  -«2)/2p2l 

d4x  ,  2 

— jr  +  x  -h/2p 

L  J 

dz 

,8 

d  x 

dz 

[4p2/<‘-2’]  £ 

=  4p2h/(l-e2) 

A  complete  solution  of  equation  (17)  is 


(17) 


3  2  2 

x(z)=A  sin<rz+Bcoso-z+Csinh(rz+Dcosh<rz+Ez  +Fz  +Gz+H  -hz  /24 

<r4=4p2/(l-€ 2)  (18) 


The  eight  constants  of  integration  are  obtained  from  the  boundary  conditions 
(3)  and  (4).  Once  again,  the  time  dependency  factors  out;  this  would  not 
have  happened  had  damping  been  present  so  that  a  solution  more  general 
than  (15)  would  have  to  have  been  used.  These  reduced  boundary  conditions 
in  x  (z)  and  y  (z)  are 
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(l-o 

d3x 

h  KjX  - 

2m  jp2 

o 

II 

I 

£ 

(l+o 

.3 

<i  y  . 

dz 

+  kiY  + 

2m  ^p2 

(x-y)  =  0 

d2x 

Z?  = 

fi 

=  0 

-(1-0 

d3x 

^7 

+  K2x  - 

2m2p2 

(x-y)  =  0 

-(l+o 

dzJ 

+  K2y 

f  2m2p‘ 

2  (x-y)  =  0 

d2x 

d7  = 

d2y 

17 

=  0 

(19) 


(20) 


When  (18)  is  put  into  (19)  there  results  eight  simultaneous  linear  nonhomo- 
geneous  equations  for  the  constants  of  integration.  For  general  values  of 

р,  «  ,K.  and  m,  the  constants  will  be  bounded  functions  of  these  system 
parameters.  However,  for  certain  combinations  of  parameters,  the  amp¬ 
litudes  will  become  unbounded;  this  happens  when  the  determinant  of  the 
coefficients  of  the  linear  equations  vanishes.  These  values  of  p  for  fixed 

с ,  m,  and  K  are  the  secondary  or  gravity-induced  critical  speeds  and  are 
given  by 

<r6  <l-t2)2  (coscrcoshtr -1)  +  <r3  (l-*2)  [jnij+m^l-c  2)  ^-(Kj+Kj) J  x 

x  jsinhacos  ir-cosh<rsino-J  -2  |mj(l-*  2)  7-KjJ  £m2(l-«  2)  7 -K2J  x 

xjsino-  sinhoj  =0  (21) 
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It  can  be  seen  from  comparison  of  Eq.  (21)  and  the  first  of  (12a) 
that  these  curves  are  identical  if  QC.  is  replaced  by  <r  and  c  is  replaced 
by  t  .  Thus  from  the  plot  of  A.,  in  the  p~c  plane,  we  can  obtain  the 

1  2 

secondary  critical  speed  by  replacing  p  with  2p  and  <  with  c  .  Thus 
the  gravity  critical  speeds  for  the  round  shaft  are  half  the  natural  frequen¬ 
cies  of  the  shaft. 
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EFFECT  OF  MASS  ECCENTRICITY 


It  is  a  practical  impossibility  to  manufacture  a  shaft  which  is  per¬ 
fectly  straight  and  uniform  throughout  its  cross  section.  Consequently, 
there  exists  some  small  mass  eccentricity  along  the  length  of  the  shaft  even 
after  it  has  been  balanced.  As  the  shaft  rotates,  the  small  centrifugal  forces 
may  be  sufficient  to  cause  vibrations  which  render  the  system  unserviceable. 
If  f  (z)  and  g  (z)  are  the  eccentricities  along  the  principal  directions,  the 
equations  of  motion  become 


(1-C) 

B4u 

.  n2  /*2u 

p  (r? 

=  P2f(*) 

(i+e) 

*4v 

,  2  .^2v 

+  P  (/rrz 

=  p2  g  (*) 

Note  that  in  the  rotating  coordinate  system,  the  centrifugal  forces 
are  constant  while  in  a  fixed  system  these  forces  are  periodic  with  a  fre¬ 
quency  equal  to  that  of  the  shaft.  A  possible  motion  of  the  shaft  is  a  "frozen 
whirl"  where  it  is  deformed  but  does  not  vibrate.  This  implies  that  there 
is  no  time  dependency  of  the  displacements  along  the  principal  directions; 
therefore  we  try  a  solution  of  the  form 

u  =  u  (z)  v  =  v  (z)  (23) 

With  the  formulation  (23)  the  equations  (22)  become  uncoupled.  The 
actual  form  of  the  eccentricities  insofar  as  stability  is  concerned  is  imma¬ 
terial.  The  solution  of  these  equations  is 

u  -  A  sina^z  +  B  cosa^z  +  C  sinha^z  +  D  cosha^z  +  ff  (z) 
v  =  E  sina^gZ  +  F  cos o'gZ  +  G  8inh*gZ  +  H  cosho^z  +  <p(z)  (24) 
*x4  =  p2/(l-E)  4*  =  p2/(l+*) 
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where  f  (s)  and  jt{ s)  are  particular  integrals  of  Eqe.  (22).  The  constants 
of  integration  are  found  from  the  boundary  conditions  (3)  and  (4).  This  re> 
suits  in  two  sets  of  four  non-homogeneous  simultaneous  linear  equations. 

The  determinants  of  coefficients  of  these  sets  is  identical  with  those  given 
by  equation  (12).  This  should  have  been  expected  since  (12)  was  developed 
on  the  basis  of  vanishing  frequency  which  is  precisely  the  "frosen  whirl"  of 
Eqs.  (23). 

It  is  apparent,  then,  that  the  flattened  shaft  possesses  two  sets  of 
natural  frequencies  like  a  non- symmetric  beam  which  vibrates  in  two  mu¬ 
tually  perpendicular  directions.  Moreover,  all  frequencies  between  these 
iiaj'acent  criticai  speeds  are  also  critical.  For  the  symmetric  shaft,  these 
natural  frequencies  will  coincide. 
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ROUND  SHAFT  IN  UNSYMMETRIC  BEARINGS 

In  most  practical  cases,  the  shaft  will  be  symmetric  while  the 
supports  possess  unequal  flexibility  in  two  mutually  perpendicular  directions. 
In  determining  the  natural  frequencies  of  the  system,  we  do  not  need  to  use 
the  previous  notation.  Since  the  equations  of  motion  and  the  boundary  con¬ 
ditions  are  uncoupled  it  is  necessary  to  consider  motion  only  in  one  direction. 
The  equation  of  motion  for  the  undamped  system  is 


while  the  boundary  conditions  are 


at  Z  s  0 


+  kj  U  +  mj 


(26a) 


at  Z  =  L  -El  t  k,  U  +  m_  =  0  (26b) 

'dz*  * 

ifu  s  0 
Vz2 

Here  El  is  the  flexural  rigidity  of  the  shaft  while  the  constants  kj  and 
are  the  stiffnesses  of  the  supports  at  Z  s  0  and  Z  =  L  respectively. 

A  solution  of  Eq.  (25)  is 


where 


U  =  (A  coshp  £  +  B  cosp  j-  +  C  sinp  £  +  D  sinp  j-  )e  (27a) 


,4  _  *X2 L4 


The  constants  of  integration  are  found  from  the  boundary  conditions  (26). 

A  non-trivial  solution  for  these  constants  results  in  the  frequency  equation 
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In  terms  of  the  redefined  non-dimensional  parameters 


Mi 

2  _  \  2  If  L  «4 

P  - x  -ET  p 

1  "  “El” 

mi  =  jr 

k2L3 

M2 

(29) 

2  3  nr- 

m2  =  7T 

the  frequency  equation  (28)  becomes 


P^  (1-cos  p  cosh  P)  +  p3  +  K2)-(mj  +  m2)  p4~||^sinh  P  cos  p  -  cosh  p  sin  P 


+  2  "  Ki  J  ^  -  K2  J  sinpsinhp  =  0  (30) 

For  the  special  case  of  identical  bearings  Kj=K2=K,  m^m^^m  and 

„  ft4  p3  cosh  p  sin  p  -  sinh  P  cos  p  i  (sin  p  -  sin  p) 

=  mp  T  . . .  y  a  '  'fl  11  . .  . 


2's'l'n  P  'sinK  P 

Taking  the  upper  and  lower  signs  separately,  we  obtain  the  frequency 
equations 


(31) 


K  —  m  p  +  (tanh  ^  t  tan  ^ ) 

3 

K  =  m  p*  +  (coth  ^  ”  cQt  \  ) 


(32) 


The  natural  frequencies  \ 


for  specified  values  of 


the  system  parameters  m  and  K  can  be  found  by  a  simple  graphical  pro¬ 
cedure.  The  functions 
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GRAPHICAL  SOLUTION  FOR  NATURAL  FREQUENCIES 


are  plotted  in  Fig.  2.  For  particular  values  of  K  and  m  we  need  only 
4  4 

plot  m tt  <t  as  indicated  and  shift  this  curve  upward  through  a  distance  K 
(dashed  curve  in  Fig.  2).  The  intersections  of  this  curve  with  the  curves 
F+  and  F"  give  those  values  of  p  which  satisfy  the  frequency  equation. 
Note  that  these  curves  asymptotically  approach  p  =  njr(n  =  1,  2,  ... )  as  K 
approaches  infinity.  Thus,  as  expected,  X  =(n^r^/L^)^I /if')  are  the  natural 
frequencies  for  a  pinned-pinned  shaft  in  rigid  self-aligning  bearings.  As  a 
further  check,  we  note  that  the  intersections  of  these  curves  with  the  0L- 
axis  (m  -  0,  K  =  0)  correspond  to  the  natural  frequencies  of  the  free-free 
shaft.  The  expressions,  of  course,  also  provide  the  natural  frequencies  of 
a  uniform  beam  under  the  same  type  of  end  constraints.  A  numerical  ex¬ 
ample  will  be  found  in  Appendix  D. 

The  natural  frequencies  for  various  mass  and  stiffness  ratios  are 
given  in  Appendix  E,  Figures  (El)  to  Ell).  This  appendix  also  contains  a 
treatment  of  six  other  types  of  bearing  configurations  which  are  of  practical 
importance. 
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UNSYMMETRIC  SHAFT  IN  IDENTICAL  BEARINGS 


The  graphical  procedure  developed  in  the  laat  section  can  be  ueed 
to  solve  the  frequency  equations  (13).  The  first  of  these  equations  is  solved 
for  K  in  terms  of  OL  to  obtain 


K 

T-6 


K 

TTi 


Ofi  <Vi 

(tanh  y  +  tan  y- ) 

0i\  Oi  1 

(coth  y  -  cot  y  ) 


(34a) 

(34b) 


TheTrequency  (^uatronn[T3b)  i»  i3entical  to  the  above  when  -  f  end  fit,  are 
replaced  by  +  £  and  respectively.  In  terms  of  the  physical  constants 
of  the  system,  the  parameters  in  Eqs.  (34)  are 


K  _  kL3  _  kL3  rJ  4  _  X2rL4  _  X2^4 

rrg  a  ST T7J  -ET7  *  *1  *  -STT^y  ET^“ 


m  = 


M 


7T 


K  vt3 

tf-£=  -sn+cT 


xVi4 

5(1+*) 


(35) 


Comparison  of  equations  (35)  and  (34)  with  (29)  and  (32),  respectively, 
indicates  that,  aside  from  subscripts  on  the  flexural  rigidities,  the  two  sets 
of  frequency  equations  are  identical. 

The  calculations  for  a  particular  rectangular  shaft  are  given  in 
Appendix  D.  The  dependence  of  critical  speeds  on  the  shaft  stiffness  aniso¬ 
tropy  €  is  shown  in  Fig.  3.  The  solid  lines  represent  the  first  three  gravity 
critical  speeds.  '  The  first  two  natural  frequencies  enclose  a  speed  range  for 
which  the  transient  solution  is  unstable.  The  effect  of  damping  on  this  criti¬ 
cal  speed  range  is  indicated  schematically  in  Fig.  (Cl). 
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EXPERIMENTATION 


The  teat  fixture  uaed  in  this  program  may  be  considered  as  consist* 
ing  of  three  parts:  the  drive  section,  which  was  simply  a  source  of  uniform 
rotational  motion;  the  experimental  section,  which  consists  of  a  very  stiff 
beam  upon  which  the  shaft  is  installed  in  spring -mounted  bearings,  and  the 
instrumentation,  by  means  of  which  the  rotational  speed  and  forces  on  the 
test  shaft  bearings  are  measured.  The  mechanical  parts  may  be  seen  in 
Fig.  4,  and  the  electrical  and  electronic  parts  in  Fig.  5. 

The  drive  system  has  as  its  principle  element  a  28  v.  d.  c.  motor 


ation  and  wide  running  speed  range.  Power  is  supplied  to  it  by  a  220  v  trans¬ 
former  seen  on  the  floor  behind  the  oscilloscope  cart  in  Fig.  5  and  by  two 
rectifiers  and  two  variacs  shown  at  the  left  of  Fig.  5.  (It  was  determined 
that  fluctuations  in  the  110  volt  line  were  too  great  to  allow  accurate  speed 
control. )  The  voltmeter  and  ammeter  on  the  front  of  the  large  rectifier 
case  are  in  the  armature  circuit.  Speed  is  regulated  by  adjustment  of  the 

larger  variac  on  top  of  the  case  at  the  front.  The  motor  is  coupled  to  a 

2 

flywheel  whose  polar  moment  of  inertia  is  59.  64  in.  -lb-sec  .  (The  polar 
moment  of  inertia  of  the  flat  shaft  is  0.  00505  in.  -lb-sec^. )  The  flywheel 
is  coupled  to  the  test  shaft  by  a  connecting  shaft  and  two  universal  joints. 

A  portion  of  the  connecting  shaft  and  one  of  the  universal  joints  may  be  seen 
in  Fig.  6.  The  drive  section  is  mounted  on  a  stiff,  cross-braced  table 
which  has  incorporated  in  it  a  vibration  damping  pad,  which  reduces  any 
vibration  transmitted  to  the  experimental  section. 

The  experimental  section  has  as  a  base  a  welded  box  beam  11-1/2 
in.  wide  by  19-1/2  in.  deep,  with  four  stiffening  baffles,  all  made  from  3/4 
in.  steel  plates.  Two  A-frames  are  bolted  to  the  beam  for  mounting  the 
springs  and  bearings.  The  stiffnesses  of  the  A-frames  in  both  horisontal 
and  vertical  directions  are  about  100  times  as  great  as  that  of  the  springs 
used  in  the  program,  thereby  insuring  that  the  springs  constitute  the  sole 
significant  flexibility  at  the  bearings.  Adjustable  bolts  through  the  A-frames 
have  cups  at  one  end  into  which  the  ends  of  the  springs  are  fitted.  Similar 
cups  on  the  bearing  blocks  accept  the  other  ends  of  the  springs,  thus  allow¬ 
ing  motion  in  either  the  horisontal  or  vertical  direction. 
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Self-aligning  anti -friction  bearings  arc  us  ad  to  minimi  sa  mom  ants 
appliad  to  tha  bearing  block  by  rotation  or  banding  of  the  shaft.  Axial  mo¬ 
tion  is  prevented  by  tha  fixity  of  the  connecting  universal  joint  on  tha  test 
shaft  and  the  held  bearing  on  the  outboard  flywheel  journal.  The  test  shaft 
is  not  held  in  either  bearing.  The  bearing  blocke  provide  mounts  for  the 
accelerometers. 

Speed  sensing  was  effected  by  means  of  a  light  source  and  an  elec¬ 
tric  eye  directed  toward  a  strip  of  aluminum  tape  on  tha  flywheel  rim  (see 
Fig.  4).  Since  the  rest  of  the  rim  was  covered  with  black  paint,  the  elec¬ 
tric  eye  emitted  a  voltage  which  rose  and  fell  once  during  each  revolution. 
The  instrument  at  the  top  of  the  rack  at  the  right  in  Fig.  5  is  a  tachometer 
which  was  used  for  rough  measurements  of  the  speed.  In  addition,  the  out¬ 
put  of  the  electric  eye  was  fed  into  the  oscilloscope  where,  by  adjustment  of 
the  sweep  time  base  and  measurement  of  the  trace  length  for  one  revolution, 
accurate  speed  measurements  were  made.  Furthermore,  the  oscilloscope 
was  used  for  vibration  phase  measurement,  through  observation  of  the  rela¬ 
tive  positions  cf  the  electric  eye  trace  and  the  spike  in  the  accelerometer 
trace. 

Bearing  forces  were  sensed  by  C.  R.  L.  accelerometers  on  the 
bearing  blocks,  and  measured  by  the  Ballantine  Model  300  voltmeter  on  the 
shelf  in  the  rack.  Other  instrumentation  in  the  rack  includes  selector 
switching  for  location  and  direction  of  accelerometers  used,  calibration 
circuitry,  and  cathode  follower  isolation  amplifiers  for  the  accelerometers. 

The  experimental  procedure  was  directed  solely  toward  the  deter¬ 
mination  of  the  critical  speeds  of  rotation.  Displacements  were  not  found 
at  this  time  since  analytic  results  were  not  available  for  comparison. 

Tests  were  conducted  on  two  shaft  systems:  a  50-in.  long  1 -in. 
diameter  round  shaft  and  a  50-in.  long  7/8  by  1-1/2-in.  rectangular  shaft. 
The  bearing  flexibilities,  limited  by  the  availability  of  commercial  coil 
springs,  were  2500  and  5300  lbs/in.  An  infinite  amount  of  stiffness,  i.  e. 
rigid  bearings,  was  obtained  by  compressing  the  springe  to  the  bottoming 
point.  The  balancing  fixtures  on  the  round  shaft  were  not  needed  since  the 
amplitudes  of  vibration  near  resonance  were  not  severe  as  long  as  the  rotor 
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passed  quickly  through  this  critical  speed.  No  attempt  was  made  to  balance 
the  flat  shaft  because  unstable  speed  ranges  and  gravity  critical  speeds  exist 
independently  of  any  mass  eccentricity  of  the  shaft. 

The  results  for  the  round  shaft  are  presented  in  Pig.  7  and  8  in 
terms  of  the  horiaontal  and  vertical  accelerations  at  the  bearing  remote 
from  the  driving  end.  The  observed  first  resonant  speed  was  at  1640,  1720, 
and  1770  rpm  for  bearings  of  2500  lb/in. ,  5300  lb/in.  and  infinite  stiffness, 
respectively.  The  secondary  gravity  critical  speeds  were  not  noticeable. 

The  resonant  speeds  were  apparent  from  the  large  vibrations  and  noise  levels. 

The  theoretical  predictions  of  Appendix  D  are  1720,  1870,  and  1910 
rpm.  The  calculations,  however,  did  not  include  the  1.  15  lb  weight  of  the 
balancing  fixtures.  The  corrected  frequencies,  based  on  a  new  mass  of 
shaft  and  elastic  modulus  E  =  29  x  10^  psi,  are  1610,  1750,  and  1790  rpm, 
which  are  in  very  good  agreement  with  the  experimental  results. 

The  results  for  the  flat  shaft  with  bearing  stiffness  of  2500  lb/in. 
are  presented  in  Fig.  9*  In  this  case  the  apparent  secondary  critical  speeds, 
occurring  at  950,  1450  rpm  were  quite  severe.  In  the  rough  running  range 
from  1600  to  2100  rpm  it  was  difficult  to  maintain  constant  speed  because  of 
the  large  vibrations  of  the  rotor.  At  approximately  2100  rpm  the  shaft  be¬ 
came  disengaged  from  the  bearing  mounts  causing  some  minor  damage  to 
the  apparatus. 

The  analysis  of  Appendix  D  predicted  gravity  critical  speeds  at  910, 
1500,  and  2140  rpm  and  a  critical  speed  range  from  1600  to  2120  rpm.  (For 
E  =  29  x  10^  psi  the  predictions  would  be  900,  1470,  and  2100  rpm  for  the 
gravity  critical  speeds  and  1570  to  2080  rpm  for  the  critical  speed  range. ) 

It  is  clear  that  the  destructive  vibrations  near  2100  rpm  were  caused  when 
the  effects  of  gravity,  mass  eccentricity,  and  shaft  stiffness  inequality  were 
felt  at  the  same  time. 
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CONCLUSIONS 


The  natural  frequencies  of  a  uniform  undamped  shaft  rotating  in 
heavy  undamped  flexible  bearings  have  been  found  for  two  specific  cases: 

1)  a  shaft  with  non-symmetric  cross  section  rotating  in  heavy,  self-aligning 
bearings  whose  principal  flexibilities  arc  equal,  and  2)  a  symmetric  shaft 
rotating  in  heavy  bearings  of  unequal  principal  flexibilities.  Secondary 
critical  speeds  induced  by  gravity  are  also  determined.  The  general  case 
of  non-symmetry  in  both  shaft  and  bearings  is  not  amenable  to  solution. 

In  the  first  case,  a  doubly  infinite  set  of  natural  frequencies  are 
found.  The  spread  between  corresponding  frequencies  of  each  set  increases 
from  zero  with  increasing  shaft  anisotropy.  Centrifugal  forces  induced  by 
mass  eccentricity  excite  a  steady  state  response  at  all  speeds  except  the 
natural  frequencies  where  the  response  becomes  unbounded.  The  transient 
solution,  however,  will  become  unbounded  at  any  speed  which  lies  between 
corresponding  natural  frequencies.  Damping  in  the  system  will  not  com¬ 
pletely  eliminate  this  range  of  unstable  operating  speeds.  The  unstable 
nature  of  this  transient  makes  it  imperative  that  this  type  of  shaft  design 
be  avoided  wherever  possible. 

For  the  case  treated  above,  the  secondary  gravity  critical  speeds 
can  be  found  with  little  extra  effort  from  the  data  on  the  natural  frequencies. 
For  slight  shaft  asymmetry  these  critical  speeds  are  approximately  one-half 
the  natural  frequencies. 

For  the  special  case  of  identical  bearings,  the  transcendental  fre¬ 
quency  equation  is  sufficiently  reduced  so  that  by  means  of  a  simple  graphi¬ 
cal  construction  natural  frequencies  may  be  found  for  any  combination  of  the 
remaining  system  parameters. 

When  the  symmetric  shaft  rotates  in  unsymmetric  bearings,  there 
are  again  two  sets  of  natural  frequencies,  but  no  ranges  of  critical  speeds. 

In  this  case,  there  is  a  direct  analogy  to  the  beam  whose  end  supports  possess 
different  horizontal  and  vertical  flexibilities;  the  natural  frequencies  of  each 
system  are  identical.  These  frequencies  can  be  found  from  the  same  graphi¬ 
cal  construction  mentioned  above. 
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It  is  interesting  to  compare  the  system  considered  in  this  report 
with  the  idealised  case  of  immovable  pinned -pinned  bearing e.  In  terms  of 

the  bearing -shaft  flexibility  ratio  K  and  maas  ratio  m  the  following  con¬ 
clusions  can  be  made:  1)  for  Urge  K  and  small  m  the  first  two  or  three 
natural  frequencies  are  approximately  equal;  2)  for  Urge  K  and  m  or 
medium  K  small  m  only  the  first  natural  frequency  corresponds  to  the 
idealised  case;  3)  for  small  K  the  frequencies  are  very  sensitive  to  the 
mass  ratio.  Moreover,  since  there  is  another  set  of  natural  frequencies 
corresponding  to  the  different  bearing  flexibility  in  the  perpendicular  direc¬ 
tion,  the  similarities  to  the  idealised  case  are  further  obscured.  If  the 
principal  bea ring  flexibiliti e a  are  nearly  equal  the  natural  frequencies  may 
occur  in  clusters. 

It  is  apparent  that,  aside  from  rotating  systems  whose  supports 
are  very  inflexible,  the  vibrations  of  flexible  rotors  operating  at  super¬ 
critical  speeds  will  be  sensitive  to  the  bearing  characteristics.  It  appears 
that  the  bearing  with  adjustable  stiffness  may  be  suitable  as  a  resonance 
changer,  thereby  effecting  a  decrease  in  noise  and  vibration  levels. 
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In  summary, 

1.  Bearing  mass  and  flexibility  have  considerable  effect  on  the  critical 
speeds  of  flexible  rotors.  Generally  only  the  first  and  second  nat¬ 
ural  frequencies  of  flexural  vibrations  are  accurately  predicted  by 
the  idealized  assumption  of  rigid  (immovable)  bearings. 

2.  When  the  bearing  stiffnesses  in  two  mutually  perpendicular  direc¬ 
tions  are  unequal,  the  occurrence  of  critical  speeds  is  doubled. 

It  is  possible  for  these  speeds  to  occur  in  clusters  thereby  giving 
the  appearance  of  a  range  of  unstable  operating  speeds. 

3.  Shafts  of  non-symmetric  cross  section  should  be  avoided  whenever 
possible  because  of  instabilities  which,  in  general,  cannot  be  re¬ 
moved  by  damping. 

4.  Half  critical  speeds  in  symmetric  shafts  are  unimportant. 

5.  Bearings  with  adjustable  stiffness  offer  a  means  of  changing  reso¬ 
nance  thereby  allowing  smooth  operation  in  between  critical  speeds. 
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RECOMMENDATIONS  FOR  FURTHER  WORK 

In  this  report  only  the  effects  of  bearing  mass  and  stiffness  on  the 
natural  frequencies  of  a  uniform  shaft  were  found.  Account  should  also  be 
taken  of  shafts  of  non-uniform  cross  section  which  may  be  carrying  thin 
heavy  disks  on  interior  and/or  overhanging  sections.  It  is  also  of  practical 
interest  to  determine  the  magnitude  of  vibrations  for  various  combinations  of 
shaft  system  parameters. 

In  addition,  the  possibility  of  using  flexible  bearings  as  a  means  of 
controlling  vibrations  should  be  pursued.  The  chief  cause  of  excessive  vi¬ 
bration  is  an  unbalanced  rotor  operating  near  a  critical  frequency.  Preci¬ 
sion  balancing  and  conservative  design  are  two  common  methods  of  vibration 
control.  While  it  is  standard  procedure  to  balance  all  critical  rotating  ele¬ 
ments,  these  rotors  may  eventually  operate  poorly  at  or  near  the  natural 
frequencies  of  the  system.  On  the  other  hand,  conservative  design  in  many 
cases  would  be  too  costly  and  impractical  especially  if  the  system  operates 
over  a  wide  range  of  speeds.  Adjustable  bearings  would  perform  the  same 
task  as  a  redesign  since  they  would  relocate  the  natural  frequencies,  i.  e.  , 
they  act  as  a  resonance  changer.  This  method  may  not  reduce  vibrations ; 
but  simply  avoid  them  by  adjusting  the  troublesome  critical  speeds  to  be 
sufficiently  removed  from  the  operating  speed. 

A  program  to  accomplish  the  above  recommendations  would  contain 
the  following  elements: 

1.  Develop  the  dynamics  of  non-uniform  shafts  carrying  thin 
heavy  disks.  This  would  include  gyroscopic  effects  and 
the  effects  of  bearing  mass,  damping,  and  stiffness. 

2.  Determine  the  amplitudes  of  vibration  for  various  shaft 
system  parameters  with  particular  emphasis  on  changes 
of  bearing  stiffness. 

3.  Perform  experiments  which  will  verify  the  analytical 
results, 

4.  Evaluate  the  feasibility  of  adjustable  bearings  as  a 
resonance  changer. 
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APPENDIX  A 


THE  EQUATIONS  OF  MOTION 


It  is  convenient  to  derive  the  equations  of  motion  in  the  rotating 
Vj,  V 2  coordinate  system  (Fig.  Al).  In  determining  the  boundary  condi¬ 
tions,  we  shall  make  use  of  the  fixed  X,  Y  coordinate  system.  We  develop 
the  equations  in  their  most  general  form  and  then  specialize  to  the  problem 
at  hand. 


Let  Sj  (Z)  and  S2  (Z)  be  the  principal  stiffnesses  in  the  and 
directions,  respectively.  Classical  beam  theory  yieldB  the  following 
relations  between  shear  force  F,  bending  moment  B  and  displacements 
Vj,  V2  of  the  neutral  axis: 


Fl  =  - 


*B2 


F2=  - 


BBj 

3T 


B2  =  Sl 


B  z 


bi=s2 


b2v 


DZ 


(A-  1) 


We  recall  that  the 
tern  which  rotates 

al  = 


the  equations  of 

2  ^V2 

-  srv.  -  zn—A) 

1  ~z> 1 

(A-3) 

2  J>V1 

.  nZV2+  2Qyr) 


For  a  horizontal  shaft  of  density  tf(z)  per  unit  length, 
motion  are 

^F 

-  Y(Z)  g  cos  n  t  +  — J-=  Y(Z)  ( - i- 

^  z  dr 

"dF, 


J  (Z)  g  Sin  Q  t  +  =  y(Z)  (  — i. 


components  of  acceleration  referred  to  a  coordinate  sys- 
at  a  constant  speed  n  are 


2vi  2  Bv2 

— J-  -  n2v  -  2n  z 

3t2  1 


I A-?.\ 


a2vz 


2  "avi 

n  V2  +  20 
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where  the  first  term  of  each  of  equations  (A-3)  represents  the  gravity  force. 
Substituting  relations  (A-l)  into  (A-3)  we  obtain 


[S1  <z>  Tjr]  +  i r<z>  \j^r  -  “S  - 20  |rj  *  -  r<z) 

[s2(z)  -Jfj  *  -°S+2°^rJ  -  -Y<z 


g  cos  nt 


(A -4) 

Y^z)  g  sin  ft  t 


When  the  shaft  is  uniform,  Sj,  S2.  and  Y  are  constants  and,  after  suitable 
non-dimensionalization,  equations  (1)  and  (2)  result. 

In  the  rotating  system,  the  component  of  force  exerted  by  the  shaft 
on  the  bearing  is 


-  S  -2—L  at  Z  =  0  and  +  S  — 


at  Z  =  L  . 


(A-5) 


The  displacements  in  the  fixed  coordinate  system  are  related  to  and 
by 


X  =  V  j  cos  ft  t  -  V2  sin  ft  t 


Y  =  V,  sin  ft  t  +  V„  cos  ft  t 

1  c. 


(A-b) 


Fx  and  F  ,  the  components  of  a  force  in  the  fixed  system,  and  F^  and  F2> 
its  components  in  the  rotating  system,  are  related  by 


F,  =  F  cos  ft  t  +  F  sin  ft  t 
lx  y 


F-  =-F  sin  ft  t  +  F  cos  ft  t 
2  x  y 


Now  the  force  exerted  by  the  springs  on  the  mass  is 


(A-7) 


F  =  -k  X;  F  =  k  Y 

xx  y  y 


(A-8) 
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where  k  ,  k  are  the  flexibilities  of  the  springs  in  the  X  and  Y  directions, 
x  y 

respectively.  Substituting  (A-6)  and  (A-8)  into  (A-7)  we  obtain 


as  the  components  of  the  force  of  the  spring  on  the  bearings. 


Finally,  when  (A-2),  (A-5),  and  (A-9)  are  taken  into  account,  the 
equations  of  motion  of  the  bearings  become 

at  Z  =  0 
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We  see  that  these  end  conditions  have  time  dependent  periodic  co¬ 
efficients  when  the  flexibilities  k  and  k  at  the  bearings  are  not  equal. 

x  y 

When  referred  to  the  fixed  coordinate  system,  the  boundary  conditions  have 
periodic  terms  arising  from  the  shaft  anisotropy.  Moreover,  the  equations 
of  motion  in  the  fixed  coordinate  system  also  have  time-dependent  periodic 


coefficients.  For  the  special  case  of  equal  bearing  flexibilities 

and  k„  =  k_  =  k_  the  nondimensionalized  equations  (3)  result. 
2x  2v  2 
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APPENDIX  B 


AN  EXAMPLE  OF  AN  UNSTABLE  TRANSIENT 
OF  A  ROTOR  CONTAINING  DAMPING 


The  transient  solution  of  the  damped  non- symmetric  shaft  systems 
considered  in  this  report  can  be  shown  to  be  unstable  for  certain  ranges  of 
operating  speed;  however,  for  the  sake  of  clarity  we  shall  demonstrate  this 
behavior  for  a  simplified  lumped  parameter  system.  Consider  an  unsym- 
metric  shaft  of  negligible  weight  supporting  a  mass  m.  In  terms  of  an  x, 
y  coordinate  system  which  rotates  with  the  shaft,  the  equations  of  motion 


•  •  ~y  •  • 

m  (y  -  D  y+2Qx)  +  C1  y  +  kpr  =  0 


m  (x  -  SI x  -  2  n  y  )  +  C^  x  +  k^x  =  0 


(B-l) 


Here  fi  is  the  constant  angular  speed  of  the  shaft,  C  is  the  damping  factor, 
and  k  is  the  shaft  stiffness. 

Since  the  differential  equations  (B-l)  have  constant  coefficients, 


the  solution  has  the  form 


y  —  A  c 


...n.H 


•  »»» 


The  frequency  equation  for  X  is  found  from  the  requirements  that  the  de¬ 
terminant  of  coefficients  of  A  and  B  vanish.  This  results  in 


0  =  mXZ  +  Cj\+  (kj  -  mn2)  2S?Xn 


(B-3) 


-2nXm 


mX2  +  C2X  +  (k2-mn2) 


m2  X4  +  m(Cj  +  C2)  X3  +  |^(kj  +  m  +  ^1^2  +  J  X2 

+  [Cl  (K2  -  mf22)  +  C2  (kj  -  m  X  +  (kj  -  m Si2)  (k2  -  mfl2)  =  0 
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In  general,  the  solution  of  the  frequency  equation  (B-3)  has  two 
sets  of  complex  conjugate  roots.  If  the  solution  (B-2)  is  bounded,  the  real 
part  of  the  complex  frequency  X  must  be  less  than  or  equal  to  zero.  When 
Real  {X)  >  0  the  transient  is  unstable  since  it  increases  exponentially  with 
time. 

While  it  is  possible  to  obtain  a  relation  among  the  shaft  system 
parameters  for  which  Real  (X)  >  0,  it  would  be  too  cumbersome  to  be  in¬ 
formative.  It  is  found  that  critical  speeds  occur  in  a  range  whose  width 
depends  on  the  stiffness  difference  (k^k^/^+kj)  and  the  damping  factor  C. 
This  critical  speed  dependence  is  indicated  schematically  in  Fig.  Bl. 

Since  the  right-hand  sides  of  Eqs.  (B-l)  are  zero,  the  casual  ob  ¬ 
server  may  conclude  that  the  system  is  in  free  vibration.  One  may  wonder 
why  the  homogeneous  solution  of  a  damped  system  can  be  unstable  when  it 
is  known  that  any  freely  vibrating  dissipative  system  will  eventually  come 
to  rest.  However,  the  equations  of  motion  (B-l)  are  really  forced  since 
the  angular  speed  is  kept  constant.  At  the  critical  speeds,  more  energy 
is  added  to  the  shaft  than  is  absorbed  by  dissipative  elements.  This  can 
happen  only  when  there  exists  a  shaft  stiffness  inequality  which  some  authors 
have  considered  apparently  equivalent  to  "negative  damping". 
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APPENDIX  C 


DERIVATION  OF  CRITICAL  SPEED  EQUATION 
FOR  A  NON-SYMMETRIC  SHAFT 


A s  |i  approaches  aero,  the  conditions  (10)  and  (11)  become 

X  Aj  [(l-tf)A'3  +  Kj-mjp2]  +  X  fAj  (l-e)rf3  ♦  KrmlP2  (jj|  )1 
j  »  1  L  J  j  -  1  L  J 

4  p  i  6  ^  i  p  “j 

X  Aje  +  X  Aje  -(l-^j3+K2-m2p2  (^|)  - 

j  =  1  L  J  j  =  1  L 


w 

X 

j  *  1 

8  <*j 

Z>j2«  a<  *  ° 

j  =  1 


(C-l) 


TTT  I  Aj  ni+O^fK^mjP*] 

j  a  5  ^ 

Tit  Y.  . 

j  «  5  L  J 


T7i  ^  °^j2  Aj  "  0 
J  *  5 


-€ 

TTc 


8  i 

^  OC  j2  e  A.  ■ 
j  -  5  * 
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where 


Mr,  l1/4 

«  (i)  j_1  |_P2/(1  -6)  J 

j  =  1,2,  3.  4 

.  _r  .  n1/4 

*J 

*  (i)  3  |_pV(1  +6)J 

j  =  5, 6,  7.8 

If  A,  the  determinant  of  coefficients  of  equations  (C-l).  were  written  out 
explicitly,  we  would  note  that  the  elements  of  the  lower  left  four  by  four 
sub-determinant  vanish.  Consequently,  the  eighth  order  determinant  is 
expressible  as  the  product  of  two  fourth  order  determinants;  that  is  (see 
page  C-3) 


256  6 
"  ‘  TTg 


The  expression  for  6^  is  obtained  from  6^  when  (1-6)  and  are 
replaced  by  (1+6-)  and  respectively.  By  a  simple  rearrangement 
of  columns  of  (C-2),  we  can  show  that  6  j  and  6^  *re  respectively  equi¬ 
valent  to  Aj  and  A^  of  equations  (12). 
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APPENDIX  D 
NUMERICAL  EXAMPLES 


I.  ROUND  SHAFT 

We  are  to  find  the  natural  frequencies  of  a  50>in.  long  1-in.  dia¬ 
meter  steel  shaft  mounted  in  self-aligning  bearings.  The  6.  2  lb  bearing 
blocks  are  connected  to  a  rigid  foundation  by  means  of  springs  of  2500  lbs/in. 

stiffness. _ The  total  mass  of  the  springs  is  negligible  in  comparison  to  the _ 

mass  of  the  bearing.  The  physical  constants  of  the  system  are 

L  =50  in.  length  of  shaft 
y  gL  =  1 1.  1 1  lb  weight  of  shaft 

Mg  =  6.  2  lb  weight  of  bearing 

k  =  2500  lb/in.  spring  stiffness 

E  =  30  x  10^  lb/in.  Z  Young's  modulus 

8  =  5. 75  x  10  lb  sec  /in.  mass  per  unit  length  of  shaft 

I  =  w  r  /4  *  4.91  x  10  in.  area  moment  of  inertia 

K  =  kL^/EI  =  212  bearing -ahaft^stiffness  ratio 
m  =6.  2/11.  11  =  0.  56  bearing-shaft  mass  ratio 

The  equations  for  the  natural  frequencies  X  are  given  by  ■ 

3  3 

K  -  m  r4  <r4  =  (tanh  ZZ  +  tan  ) 


where 


3  3 

K  -  m  r4  V4  =  I_JL_  (coth  2  -  cot  2  ) 


(D-l) 


v  =  h 

w 


“ET 
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In  Fig.  D-l  the  functions  given  by  Eqs.  (D-l)  and  the  dashed  curve 
212  -  0.  56  ir*  <r*  are  drawn.  The  intersections  of  these  curves,  points 
A,  B,  C,  and  D,  represent  the  first  four  dimensionless  natural  frequencies 
<r  =  0.  95,  1.31,  1.43,  and  2.  09,  respectively.  The  frequencies  X  (in 
radians  per  second)  are  given  by 


thus 

Xjl  =  180  rad/sec  =  1720  rpm 

=  342  rad/sec  =  3270  rpm 

Xj  =  409  rad/sec  =  3900  rpm 

X4  =  874  rad/sec  =  8350  rpm 
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II.  RECTANGULAR  SHAFT 


We  determine  the  rough  running  speeds  of  a  rectangular  shaft 
mounted  in  the  flexible  bearings  of  the  previous  case.  The  constants  of 
the  system  are 


when 


cross  section  a  0.  875  in.  x  1.8  in. 

r  gL  a  18. 6  lb 

Ij  a  0.  0837  in.  4 

I2  a  0. 246  in. 4 

SaEdj  +  i2)/2  «  4.45  X  106  It.  in 

6  =  (I  -I  )/d  +1  )  *  0.49 

K  kL^/S  =  61.3 

K/(l-«)  a  124 

K/(l+<)  =  42.3 

K/(1-«S  a  83.2 

m  a  0. 333 

The  frequency  equation  for  an  anisotropic  shaft  is  given  by 


K  -  mo£4  a  (tanh  +  tan  9^. ) 

k'  -  mt*4  a  (coth  -|=  -  cot  ) 


(D-2) 


w'  _  K  4  _  X2*L4  _  XZtfL4 

K  "  TTg  »  <X-  ~  ■BTTTT  "  EI7" 


when 


K  ^4  _  X2^4  _  xVl4 

it?  ’  06  -5rn?T  "Tn 
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A  graphical  solution  of  these  equations  is  obtained  by  constructing  the 
functions 


f+  (y  =  llgt  (tanh  Zgr  +  tan 
3  3 

f-  (^)  =  *-* c.  (coth  -  cot 

and 


. K*  -  m  . - -  -  -  - 

i 

The  intersections  of  the  dashed  curves  K  *  124  and  42.  3  with  the  full 
curves  in  Fig.  (D-2)  determine  the  dimensionless  natural  frequencies 


and 


respectively. 


Taking  k'  a  124  we  obtain 


a  * 

0.91 

or 

X  = 

167  rad/sec  ■ 

1600 

rev/min 

t  • 

1.25 

or 

X  = 

315  rad/sec  = 

3010 

rev/min 

i s 

1.43 

or 

X  * 

413  rad/sec  = 

3940 

rev/min 

For  K  *  42. 3  we  find 


*1  * 

0.79 

or  X 

=  222  rad/sec  * 

2120 

rev/min 

*1  a 

0.  95 

or  X 

=  312  rad/sec  * 

2980 

rev/min 

1  * 

1.27 

or  X 

s  558  rad/sec  * 

5330 

rev/min 
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In  obtaining  the  gravity  critical  speeds,  wa  can  aaaily  a  how  that  tha 
frequency  equation  is  identical  with  Eq.  (D-2)  where 

K'»K/|1*£2)  and  H  ~ 

S(l-f) 

Thus  from  Fig.  (D-2)  with  K'  =  83.  2  and 


we  obtain 


0.88  or  X  >  95  rad/aec  a  910  rev/min 

=  1.  13  or  X  =  156  rad/aec  a  1500  rev/min 

=  1. 35  or  X  *  224  r/td/aec  a  2140  rev/min 

This  particular  shaft  system,  then,  has  the  following  characteris¬ 

tics:  discrete  rough  running  speederat  910,  1500,  2140  rpm,  etc.  and  criti¬ 
cal  speed  ranges  1600-2120  rpm,  2980-3010  rpm,  3940-5330  rpm . 
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MODE  SHAPES  FOR  CASE  I 


The  three  lowest  modes  of  vibration  corresponding  to  the  pinned* 
pinned  shaft  of  case  I  are 

Uj=0.  5  cosh  2.  970  z  +  0.  5  cos  2.  970  a  -  0.  449  sinh  2.  970  s  +  5. 964  sin  2.  970  s 

u^sO.  5  cosh  4.  102  z  +  0.  5  cos  4.  102  z  -  0.  518  sinh  4.  102  z  +  0.  260  sin  4.  102  z 

UjsO.  5  cosh  4.  478  z  +  0.  5  cos  4.  478  z  -  0. 488  sinh  4.  478  z  -  0.  663  sin  4.  478  s 


These  modes  are  constructed  in  Fig.  D-3  where  the  bearing  displacement  is 
taken  as  unity.  Although  the  first  and  third  modes  are  orthogonal  to  the 
second,  they  are  not  orthogonal  to  one  another. 
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FIRST 


^  THIRD 

0 

Fig.  D3  MODES  OF  VIBRATION  K  =  212,  m  =  .  56 

D 

0 

B 


APPENDIX  E 


0 
0 

0 

f  THE  FREQUENCY  EQUATIONS  FOR  VARIOUS 

BEARING  CONFIGURATIONS 


The  frequency  equation*  for  uniform  round  shaft  system*  can  be 
obtained  in  reasonably  simple  forms.  The  various  bearing  configurations 
can  be  characterised  as  follows: 

I  Both  Bearings  Flexible  (K^  =  K^,  =  M2) 

a)  Pinned- Pinned  (already  treated) 

b)  Clamped-  Clamped 

c)  Clamped- Pinned 

II  Left  Bearing  Flexible.  Right  Bearing  Rigid 

a)  Pinned- Pinned 

b)  Pinned-  Clamped 

c)  Clamped- Pinned 

d)  Clamped-  Clamped 

The  natural  frequencies  for  a  free  type  of  end  constraint  may  be 
found  as  the  limiting  case  K-^-0,  while  the  immovable  constraint  is  found 
from  the  limit  K-^oo 
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Cate  lb  Clamped- Clamped 

The  boundary  conditions  (26)  now  become 


at  Z  =  0 


+  kjU  +  Mt 


0 


du 

VZ  = 


0 


at  Z  =  L 


M  ®*U 

^1? 


(El) 

(E2) 

(E3! 


8U 

*2  = 


0 


(E4 


Using  expression  (27)  as  the  solution  at  the  equation  o£  motion  (25),  we 
obtain  the  frequency  equation  (E6)  given  on  the  next  page.  For  the  case  o 
identical  bearings  the  frequency  equation  takes  on  the  simple  form 


K  =  mir^od^  +  2ir^o6^/(tanh  -  tan  ^£-) 

K  =  mw^oc^  +  2ir^od^/(coth  —  +  cotU<y- )  (E5 


The  first  four  natural  frequencies  for  various  mass  and  stiffness  ratios  a 
given  in  Figures  (E12)  to  (E17).  The  asymptotes  correspond  to  the  natun 
frequencies  of  immovable  bearings  which  are  found  from 


tanh 


TTOC 

~r~ 


w  °t 
~1 
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Here  the  boundary  conditions  (El),  (E2)  and  (  E3)  are  unchanged 


while  the  remaining  condition  becomes 
at  Z  ■  L  -  —  >0 


(E7) 


For  identical  bearings,  the  frequency  equation  (E9)  given  on  the  next  page 
reduces  to 

K  =  mrr*oc+  coshw^cos  »i<)+V(l-cos  ngccoah  tt«0^+4(cos  ir<<  +coshiix)2 

Ztslnh  lr  ^  c  o  i  f#C  -  ^SaOTWiatBrWfttT - 


The  first  few  natural  frequencies  for  various  mass  and  stiffness 
ratios  are  given  in  Figures  (E18)  to  (E23).  The  asymptotes  correspond  to 
the  natural  frequencies  of  the  ideal  clamped- pinned  shaft  which  are  found 
from 

tanh  t^a  tan  HOC 
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Case  Ila  Pinned  Flexible-Pinned  Rigid 


Any  number  of  frequency  equations  may  be  presented  by  varying 
the  shaft  system  parameters  K^,  K^,  and  M^.  However,  we  shall 
restrict  ourselves  to  the  case  where  one  bearing  is  practically  rigid 
(at  Z  =  L)  while  the  other  bearing  remains  flexible  (at  Z  a  0).  Accordingly 
the  boundary  conditions  (26a)  remain  unchanged  while  the  boundary  conditions 
at  the  other  end  become 


at 


Z  =  L  U  =  0, 


0 


(E10) 


where 


Again  assuming  the  solution  (27), 
4  4  ^ 

K  =  mir  +  — 7- —  ( coth  ir  oC 


1/4 


we  obtain  the  frequency  equation 
cot  ir<*  )  (Ell) 


The  first  few  natural  frequencies  for  various  K  and  m  are  given  in  Figures 
(E24)  to  (E29).  Note  that  the  idealized  case  of  pinned-pinned  immovable 
vearings  (K->oo)  furnishes  the  natural  frequencies 
<*  =  1,  2,  3,  ... 
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Case  lib  Pinned  Flexible-  Clamped  Rigid 

In  addition  to  the  boundary  condition*  (26a)  we  have 


at 


Z  =  L  U  =  0, 


8U 

TZ 


=  0 


(E12) 


The  frequency  equation  corresponding  to  the  solution  (27)  is 

K  -  mir4^4  +  ff3pd3(l  +  co*  ff  adcosh  ind) _ 

sinh  iro^cos  irod  -  cosh  ^-c^sin  IreS 


(El  3) 


2  nl/4 


The  dependence  of  the  natural  frequencies  X  on  K  and  m  is  illustrated  in 
Figures  (E30)  to  (E35). 
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Cue  lie  Clamped  Flexible -Pinned  Rigid 

The  conditions  for  a  clamped  flexible  bearing  are  given  by  conditions 
(El),  (E2),  while  a  pinned- rigid  bearing  has  the  conditions 

at  Z  =  L  U  =  0,  i-iL  =  0  (E14) 

8Z2 

The  frequency  equation  is 

K  =  m*  +  2*^ot?/(tanh  neC.  -  tan  ir<*)  (El 5) 

The  first  few  solutions  of  this  equation  are  given  graphically  in  Figures 
(E36)  to  (E41). 
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Case  lid  Clamped  Flexible- Clamped  Rigid 


at 


In  addition  to  conditions  (Ei),  (E2),  we  have 

8U 


Z  =  L  U  =  0, 


S7. 


=  0 


(E16) 


The  frequency  equation  is 


3^3, 


T,  4.4.  ir  oC  (sin  ir  o^cosh  ircc  +  cos  ir  Odsinh  iro<  ) 

K  =  mir  oC  +  - (cos  it '^cosK TrSC  -  i) -  (E17) 


The  first  few  solutions  of  this  equation  are  given  graphically  in  Figures 
(E42)  to  (E47). 


ARMOUR  RESEARCH  FOUNDATION  OF  ILLINOIS  INSTITUTE  OF  TECHNOLOGY 


Report  No.  K27 4 


E-9 


2 


3 

m  =  0.9 


5 


6 


4. 


Fig.  E1°  NATURAL  FREQUENCIES  -  CASE  la 


E-22 


t 


mii 


E-37 


8 

0 

0 

L! 

I 

L 

L 

U 

II 
11 


ssLi.-: 


D 

0 

0 

L! 

[ 

I 

L 

li 


U 

l! 

0 

li 

11 

I! 

11 

11 

II 

0 

D 


m  s  0 


Fig.  E36  NATURAL  FREQUENCIES  -  CASE  He 


E-45 


E-46 


I 


E-56 


*lr 


APPENDIX  F 


DETAIL  DRAWINGS  OF  TEXT  FIXTURE 


FI 


j  APPENDIX  F 

i 

DETAIL  DRAWINGS  OF  TEXT  FIXTURE 


FI 


\\WA'N\  A\\\\\\\\V\\V\A\'  ,v< 


a  narcn  ■ 


®®  BEARING  HOUSINS  COVER  PLATE 
nATu- Miue  »ticl  i  Ttca.o 


■pnyuk'i 


i 


